Nesting Induced Precursor Effects: a Renormalization Group Approach 



O 

o 
o 

(N 



00 
(N 



F. Vistulo de Abreu a and Benoit Dougot 6 

a Departamento de Fisica, Universidade de Aveiro, 3810 Aveiro, Portugal. 
b Laboratoire de Physique Theorique et Hautes Energies, Jussieu, 75252 Paris Cedex 05 

(February 1, 2008) 

We develop a controlled weak coupling renormalization group (RG) approach to itinerant electrons. 
Within this formalism we rederive the phase diagram for two-dimensional (2D) non-nested systems. 
Then we study how nesting modifies this phase diagram. We show that competition between p-p 
and p-h channels, leads to the manifestation of unstable precursor fixed points in the RG flow. This 
effect should be experimentally measurable, and may be relevant for an explanation of pseudogaps 
in the high temperature superconductors (HTC), as a crossover phenomenon. 

PACS numbers: 74.20.-z , 74.20.Mn, 75.10.-b 
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An important clue to understand the physics of the 
HTC may come from the unconventional normal state 
properties preceding superconducting ordering: a pseu- 
dogap is opened, which manifests through an apparent 
continuous reduction in the density of statesEJ at a tem- 
perature T*, higher than the superconducting critical 
temperature, T c . The precise definition of this inc 
ent ordering mechanism is still subject of discussion! 
Two main pictures have been proposed depending on 
whether the correlations heing built up in the pseudo- 
gap are of superconductingu or (spin, or charge) density 
wave natured. However no experimental conclusive proof 
has been given so far favoring any of these pictures. 

In the former, it is assumed that pre-formcd pairs ap- 
pear at the pseudogap energy scale T* , but coherence 
would only be established at T c . Some support in favor of 
this scenario is provided for instance by ARPES measure- 
ments. They indicate that the pseudogap shares some of 
the superconducting properties like the gap symmetry. 
Furthermore, they reveal the presence of a non- vanishing 
gap in the vicinity of the (tt, 0) point already above TjB. 
Nevertheless, it has been argued that NMR and heat 
capacity data show, that superconductivity and pseudo- 
gap are competing!], like antagonistic phases. Also, T* 
merges with T c at a critical doping in the overdoped su- 
perconducting region, when the condensation energy, the 
critical currents and the superfluid density reach sharp 
maxima. This has been argued to be irreconcilable with 
a preformed pairing scenario!. 

The amplitude of incommensurate peaks observed in 
inelastic neutron scattering varies upon doping in a way 
highly related torthe pseudogap and falls to zero in the 
overdoped regionn. ARPES has clearly shown the exis- 
tence of a Fermi surface which is lying not too far from 
the diagonals of the first Brillouin zoneQ, therefore pro- 
viding an approximate nesting. 

In the ideal case of perfect nesting, it is well known 
that density wave and-,superconducting instabilities are 
strongly interfering^ til. In this work we present a simple 
one-loop RG approach to electrons on 2D Fermi surfaces 
with flat nested regions, away from half-filling. We show 
that the strong competition arising from particle-particle 
and particle-hole corrections when nesting of the Fermi 
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FIG. 1. Left: Tomographic parameterisation of a Fermi 
surface with 3 patches per side (N c = 3): 2 on nested re- 
gions (N=2) and 1 in curved portions (M=l). The Ci (1, 1) 
is illustrated. Right: Diagrams considered. 



surface (FS) is included, naturally induces instabilities 
that can be preceded by other incipient instabilities. In 
the RG sense this appears as a crossover phenomenon: 
the RG flow passes near an unstable fixed point, captur- 
ing the physics of this fixed point at intermediate ener- 
gies, before the system orders in the final superconduct- 
ing state. This scenario naturally provides two energy 
scales and therefore bears striking similarities with the 
pseudogap phase of the HTC's. The condition of perfect 
nesting may be relaxed provided interactions are large 
enough. 

We parametrize a 2D Fermi surface with a set of 4iV c 
patches, as illustrated in figure 1. The kinetic energy 
dispersion is linearized around the Fermi surface, each 
patch being assigned a Fermi velocity, vj. In this work 
all vj are assumed equal as the results do not depend 
significantly on smooth Fermi velocity modulations. The 
most general spin rotational invariant interaction term 
which remains in an effective low energy theory only in- 
volves incoming electrons on opposite sides of the Fermi 
surfacell3't3. This remark enables us to separate the set 
of patches in two complementary subsets called right (R) 
and left (L) patches, for convenience. We may therefore 
write: 
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and 



where charge Jg(I) = J2 T ,r> '■ R i+s,t 5 t,t> Ri,t 
spin currents Jf (I) = J2t,t> '■ R t+s.r a T,r' Riy have 
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been used, and momentum variables within each patch 
have been omitted to simplify the notation. Here chain 
indices could have, been written in bold as they could 
represent vectora!3. In that case we could, for instance, 
consider a bilayer of 2D electron systems: the present 
formalism is in fact quite general allowing the treatment 
of many different systems. Eq. (Q) allows an extremely 
simple derivation of one-loop equations: only the two di- 
agrams in fig. 1, need to be considered. Requiring the 
invariance of the two-particle vertex under cutoff reduc- 
tion we arrive at (d = d/dlnD and C, = Ci/2nN c ) : 

dC c s (I, K)=J2 K, S, a) (C c a (I, K)C c s _ a (I + a,K + a) 

a, 

+3C* (I, K)C$_ a (I + a,K + a))/ (v I+a + v K+a ) 

- S ZS (I, K, 5, a)(C c a (I, K-a + 8)C<g_ a {I + a, K) (2) 

a 

+3C s a (I, K-a + 6) C' s _ a {I + a, K))/ (v I+a + v K - a+ s) 
for charge couplings, whereas for spin couplings we get: 
dCI(I, K) = Y, 5 C (I, K, S, a) (C s a (I, K)C c s _ a (I + a,K-a) 

a. 

+C c a (I.K)CI_ a (I + a,K-a) (3) 
-2CJJ, K) C s s _ a (I + a,K- a)) I (v I+a + v K+a ) - 

- S ZS (I, K, S, a)(C c a (I, K-a + 6)C s s _ a {I + a, K) 

a 

+C s a (I : K-a + S)Q_ a (I + a, K) 
+2C s a (I, K-a + S) C$_ a (I + a, K))/ (v I+a + v K - a +g) 

Here 5 C and S zs implement phase space restrictions: they 
are equal to 1 iff all electrons involved in the scattering 
are kept on the FS along the process. This happens if 
the electrons are on nested regions of the FS, or, in the 
case of the Cooper channel, if / = K. Otherwise 8° and 
S zs are set equal to zero. 

In order to characterize the possible instabilities in the 
system, we introduce several response functions, of the 
form: TZg (I,K) = -ifdte tuJt (T {O s (I;*) Oj (K;0)}) . 
For singlet superconductivity we denote TZg (I, K) = 
Sg (I, K) and the order parameter is: Ossc,S (I]t) — 
7/2 J2t t Ri+s,t (t) £j,-t (t) ■ For triplet superconduc- 
tivity, K S (I,K) = Sf(I,K), and TS c,s(I;t) = 
Y^, T Ri+s,t (i) Li t - T (t) . As interactions do not con- 
nect pairs of particles— with different S (because of mo- 
mentum conservation)!!^, the problem of finding the most 
divergent instabilities is reduced to that of diagonalizing 
matrices Sg (I, K) , with 5 fixed. The eigenvectors deter- 
mine the order parameter characterizing the instability. 

For charge density wave we define IZs (I, K) = 
Dg(I,K), with the order parameter Ocdw,s (I) = 

75 R I+S,r L hr = 75 R ll +Q ,t L M -Q ,r (i.e., 

M = /+! and Q =§)• The last definition stresses that 
interactions connect only order parameters with equal 
particle-hole center of mass, M. Thus, the problem of di- 
agonalizing Dg (I, K) is once again simplified as there is a 



fixed parameter. Finally for spin density waves we define 
similarly S dw,s (I) = 75 E r tR m+q,t l m~q,t- 

With these definitions it is straightforward n t£L ob- 
tain RG equations for the response functionsEl'Ea, of 
the form d\nH 5 {I,K) /dhxD = v with K S {I,K)= 
7r (vi+s + vi) dlZg (I, K) /dlnw. More precisely: 
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Cf_ K (K, K + S)~ 3Cf_ K (K, K + S) 

{vk+8 + Vk) 

Cj_ K (K, K + S)+ C\_ K (K, K + S) 

(v K +s + v K ) 
Cf_ K {K, I + S)+ 3Cf_ K {K, 1 + 5) 

(v 2 i-K+s + vk) 
Cj_ K {K,I + 5)-Cf_ K (K,I + S) 

(V2I-K+8 + vk) 



(4) 
(5) 
(6) 
(7) 



In this work we are concerned with the ordering of 
the electronic system. As a general trend, the effec- 
tive couplings in the one-loop approximation diverge for 
a finite value of the reduced cut-off. This is naturally 
interpreted as the signature of a low temperature in- 
stability towards an ordered state. However this di- 
vergency of the effective couplings also signals a break- 
down of the one-loop approximation. This difficulty 
may be overcome by noting that in general ratios of 
any two coupling constants reach finite asymptotic val- 
ues as the instability is approached!! In this work we 
have thus emphasized the study of these asymptotic flows 
in the space of directions for the effective coupling vec- 
tor. Representing the Cg(I,J) by a generic Cj, this 
can be done by finding the flow of non-diverging nor- 
malized couplings a = Ci/Af, with Af ~^Jj2j (Q) 2 ( tne 
sum is over all the couplings). From the RG equations 
(||) and (||) having the form dCi/dhxD = AijkCjCk 
we obtain the evolution for the normalized couplings: 
dci/du = ^2j ki Ajki Cfc q (Sij — CiCj) . We have introduced 
a new scale parameter u related to D by du = —Af 
d\nD. In this parametrization the magnitude M of 
the coupling vector diverges only as u — > 00. From the 
flow as a function of u, we can recover the more physi- 
cal variable D since dhiAf / du= Aiki Ci Ck Q which yields 
Af (it) and therefore D{u). A similar procedure can be 
applied to the computation of response functions, that 
we generically represent by TZi. Their RG equations 
have the form d In TZi /dhiD = Bik Ck , which becomes 
d In TZi /du = Bi^Ck- As u goes to infinity the r.h.s. 
reaches a finite value. This is easily seen to imply a di- 
vergency of the form TZi ~ [^ n (D / D*)]~ ai , where D* is 
the scale at which the instability occurs. The exponents 
cti are completely determined by the fixed point values 
c* of the normalized couplings. 

This procedure provides a controlled framework to 
study the complete RG flow (not restricted to the one- 
loop approximation) provided the bare couplings are 
small enough. Certainly we expect the one-loop approach 
to break down if Af is of the order unity. For a given 
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choice of initial couplings with norm A/"o, this yields a 
maximal value of u (u max ) beyond which our approx- 
imation is no longer reliable. The main point is that 
""max can be made arbitrarily large provided Mo is small 
enough. In general we expect M max ~ hi(l/-A/"o). In some 
applications we are also interested in a finite an not in- 
finitesimal A/"o regime. Higher order corrections destroy 
the fixed directions found in the one-loop approximation. 
However, we conjecture those become fixed points in the 
usual sense and that the topology of the one-loop flow 
pattern is preserved. 

As a first example of application of the present formal- 
ism we consider a non-nested Fermi surface. The spin- 
less caseJjs thoroughly discussed in the review paper by 
Shankarli3. Phase space restrictions set in through 
the 5 functions, imply that only couplings of the form 
Cs{I,I) are renormalized, due to particle-particle cor- 
rections (BCS channel). Divergence of the flow is then 
related to superconductivity. Indeed, inspired from (|J- 
|5|) , the RG flow equations can be quite simply written in 
terms of singlet Cf(I) = C%{I,I) - 3C$(I,I) and triplet 
couplings CJ(I) = C$(I, I) + CI (I, I): 

dCf(I) = a (Cj(I)) 2 dCg(I) = a (Ci(I)) 2 (8) 

where a is a Fermi velocity dependent constant. The 
phase diagram is shown in Figure 2, for a system with 
initial couplings C£(I,I) = C c and Cf(7,I) = C s . The 
only region where the flow does not diverge corresponds 
to a Fermi Liquid phasetj. In the other regions, the RG 
couplings diverge at the energy scale that we associate 
with the superconducting critical temperature, T c . 

This picture changes when we introduce nesting. As 
already discussedEjO flat regions dramatically enhance 
the number of low energy couplings to be considered. 
The occurrence of log singularities in both the particle- 
particle and the particle-hole channels induce an intricate 
interference mechanism, leading to the so-called parquet 
regime. In Figure 3, we show the new phase diagram, 
together with typical flow patterns (using the u vari- 
able, and normalized couplings) for each of the different 



phases. These results have been obtained for a partially 
nested Fermi surface with N = 6 patches along each flat 
segment and M = along each curved arc. Note that 
increasing the number of patches doesn't change signifi- 
cantly the phase-diagram. 

The modifications due to nesting are most obvious in 
the initial coupling domain corresponding to the Fermi 
liquid regime for the circular Fermi surface. For positive 
bare spin coupling, we obtain a CDW phase, as a mean 
field analysis predicts. The most striking result appears 
for repulsive couplings in the charge sector and moder- 
ate negative spin couplings. With this parametrization, 
the simple repulsive Hubbard model falls in this region. 
The flow pattern from Figure 3 clearly shows a two-step 
evolution as u is increased. Before entering the final 
asymptotic regime, the system spends a large u-time in 
the proximity of an unstable fixed point (for the nor- 
malized couplings). In this intermediate phase, the most 
diverging response is the SDW susceptibility. The cor- 
responding order parameter has the approximate form 
RpfLi + R^ 1 _ 1 L2 + ... + h.c so the most favorable nest- 
ing vector is perpendicular to the flat regions. Couplings 
involving incoming and outgoing electrons close to the 
endpoints of flat regions are the most enhanced. This in- 
termediate fixed point corresponds very likely to the one 
found by Zheleznyak et a£3 who worked with the orig- 
inal un-normalized couplings. Our approach yields the 
remarkable result that the system is finally attracted to- 
wards a stable low energy d-wave superconductor (SSCd) 
fixed point. Again, the largest couplings involve the end- 
points of the flat regions, and from the superconducting 
response function, the gap exhibits a strong anisotropy. 
It vanishes at the center of flat regions and reaches its 
maximal values at their extremities. The dramatic role 
played by these extremal regions suggests some similarity 
with the Van-Hove singularity scenario studied by sev- 
eral groupsaO'EEl. However, we haven't introduced any 
singularity in the Fermi velocity. Our model also ignores 
Umklapp processes, which are expected to play a crucial 
role in stabilizing a commensurate SDW (Neel state) , and 
in describing the proximity to a Mott insulatoillil. 
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FIG. 4. Flow of the dominant SDW (thin) and SSCd 
(thick) response functions for systems with N+M=14, but dif- 
ferent nesting qualities (N=14,12,8,3). Better nested systems 
develop instabilities at higher energy scales. SSCd response 
functions become the most diverging responses at the final 
energy scale, associated with T c . 



It is important to specify the energy scales associated 
to both fixed points. Assuming the bare bandwidth is Dq, 
the correspondence between the running scale D and the 
fictitious time u takes the form: ln(D /D) — h{u)/Mo, 
where A/"o is the length of the initial coupling vector, and 
h(u) is a function depending on the direction of this vec- 
tor. In general, h appears to be an increasing function of 
u which reaches a finite limit h(oo) as u becomes infinite. 
This defines the instability scale D c below which the sys- 
tem is superconducting. Similarly, we define a scale D\ 
for which the distance of the normalized coupling vec- 
tor to the intermediate fixed point is minimal. Let us 
note by u\ the corresponding value of u. We obtain: 
Dx/Dq = (D c /Do) h(ui)/hioD) , which holds for a chosen 
direction of the initial coupling vector. D c and D\ are 
considered here to be functions of A/"o. In practice, we 
have found that h(ui)/h(oo) is very close to unity. As 
an order of magnitude, if D c ~ 100K and Do ~ 5Q0QK, 
Di — D c ~ IK. However, SDW correlations begin to 
build up at a much higher energy scale, as can be seen 
in Figure 4. We may for instance evaluate the scale D* 
at which the SDW exponent reaches 99 percent of its 
maximal value. 
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Table 1: System with N+M=14 and normalized couplings. 

From the table above, we conclude that there can be a 
large temperature range for which a SDW enhancement 
is observed above the superconducting instability. Both 
temperature scales are seen to decrease as the length of 
nested segments is reduced. This reduction is accompa- 
nied by an increase in the minimal distance to the un- 
stable fixed point along the flow. This leads to a less 
pronounced intermediate regime, smaller exponents for 
the SDW response, and consequently a smaller magni- 
tude of SDW correlations. 

In summary, we have shown that precursor effects are 
predicted by a controlled weak coupling RG approach 
to fermionic systems. The flow approaches an unstable 



fixed point for systems with repulsive interactions, be- 
fore being attracted to the final stable superconducting 
fixed point controlling the low temperature physics of 
the system. SDW response functions are enhanced and 
dominate near the intermediate fixed point, before the 
d-wave superconducting response finally diverges. We 
believe that Umklapp processes are likely to strengthen 
this picture. Close to half-filling, they are not expected to 
compete with SDW ordering, but with superconductiv- 
ity, thereby reducing T c without affecting the cross-over 
energy scale. In spite of the weak coupling limit chosen 
here, we observe striking similarities with the pseudogap 
behavior of the HTC. This suggests that the strength of 
interactions in real materials may compensate the devia- 
tions of the observed Fermi surface from perfect nesting. 
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